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Abstract 

The resummation of logarithmically-enhanced terms to all perturbative orders 
is a prerequisite for many studies of QCD final-states. Until now such resumma- 
tions have always been performed by hand, for a single observable at a time. In 
this letter we present a general 'master' resummation formula (and applicabil- 
ity conditions), suitable for a large class of observables. This makes it possible 
for next-to-leading logarithmic resummations to be carried out automatically 
given only a computer routine for the observable. To illustrate the method we 
present the first next-to-leading logarithmic resummed prediction for an event 
shape in hadronic dijet production. 
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1 Introduction 



QCD is unique among the theories of the standard model in that both strong and weak 
coupling regimes are relevant to modern collider experiments. This manifests itself most 
dramatically in hadronic final states of high-energy collisions, whose branching pattern is 
sensitive to physics spanning the whole range of scales from the (perturbative) hard collision 
virtuality down to (non-perturbative) hadronic masses. Accordingly final states are a 
privileged laboratory for QCD studies: perturbative investigations have for example led to 
many measurements of the strong coupling, a s [1], and to tests of the underlying SU(3) 
group structure of the theory [2]; and final-states are also proving to be a rich source of 
information on the poorly understood relation between perturbative, partonic predictions 
and the non-perturbative, hadronic degrees of freedom observed in practice [3,4]. 

Among the most widely studied final-state properties are measures (v) of the extent to 
which the geometric properties of an event's energy-momentum flow differ from that of a 
Born event (the lowest order contribution to the given process, for example e + e~ — > qq). 
Fixed-order perturbative calculations, which involve a small number of additional partons, 
are suitable for describing large departures from the Born-event energy flow pattern, in 
which the extra partons are energetic and at large angles. Such configurations are however 
rare, their likelihood being suppressed by powers of the perturbative coupling. 

The most common events are instead those in which the departure from the Born 
energy-flow pattern is small, v <C 1, with any extra partons being soft and/or collinear 
to the original Born-event partons. This poses a problem for fixed-order studies because 
each power of the coupling is then accompanied by up to two powers of the large logarithm 
lnl/u, associated with soft and collinear divergences. As a result, the perturbative series 
involves terms (a s ln 2 l/v) n , and must be resummed to all orders. 

Today's state of the art calculations exploit the fact that for many measures ('observ- 
ables'), the dominant all-orders perturbative contribution can be written as an exponential 
of leading-logarithmic (LL) terms a^ln n+1 1/v. Furthermore the next-to-leading logarith- 
mic (NLL) terms, a™ ln n 1/v, factorise and can be calculated to all orders [5]. But to obtain 
this NLL accuracy one needs a detailed understanding of the observable's analytical prop- 
erties and of the corresponding phase-space integrals. Thus it is usual for an entire paper 
to be dedicated to the resummation, in a single process, of just one or two observables. 

In this letter we instead adopt the novel approach of simultaneously examining a whole 
class of observables, for which it will be possible to carry out a common analysis. The 
results, involving a 'master' formula with applicability conditions, will be relevant to a 
range of processes including e + e~ to 2 or 3 jets, DIS to 1 or 2 jets, Drell-Yan (or 7, W^, 
Higgs,. . . ) plus a jet, and hadronic dijet production. The final answer for some specific 
observable will be expressed in terms of straightforwardly (and automatically) identifiable 
characteristics of the observable. 
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2 Master formula and applicability conditions 



Let us start by taking a Born event consisting of n hard partons or 'legs' (rij of which 
are incoming), with momenta pi, ■ ■ ■ ,p n - We shall consider the resummation, in the n-jet 
limit, of n-jet infrared and collinear (IRC) safe observables — these measure the extent to 
which an event's energy flow departs from that of an n-parton event. For the resummation 
approach to be valid the observable (a function V of the final-state momenta) should: 

1. vanish smoothly as a single extra (n+l) th parton of momentum k is made asymp- 
totically soft and collinear to leg £, the functional dependence being of the form: 

(k \ at 

Vm,k)=dil^J e-<^(0). (1) 

Here Q is a hard scale of the problem; {p} represents the Born (hard) momenta after 
recoil from the emission, which is defined in terms of its transverse momentum kt 
and rapidity r\ with respect to leg £, and where relevant, by an azimuthal angle 
relative to a Born event plane. By requiring the functional form JT} (in practice, 
almost always valid), the problem of analysing the observable reduces in part to 
identifying, for each leg £, the coefficients a^, b^, d,£ as well as the function ge{4>) 
parameterising the azimuthal dependence (the normalisation may be fixed by the 
condition gi(n/2) = 1). IRC safety implies a>i > and b(, > —an (see also [6]). We 
further require the observable to be positive definite. 

2. be recursively IRC (rlRC) safe: meaning that, given an ensemble of arbitrarily soft 
and collinear emissions, the addition of a relatively much softer or more collinear 
emission should not significantly alter the value of the observable. The formal require- 
ment can be formulated as follows: we introduce momenta ftj(Aj) that are functions 
of parameters A,; such that, 

V({p},K i (X i )) = X i7 (2) 

with the condition that in the soft and/or collinear limits, Aj — > 0, the azimuthal 
angle (pi of Ki(Xi) should be fixed. Each of the momentum functions Ki(A), ^(A), 
etc. may be different as long as they all satisfy eq. (0). The conditions for rlRC 
safety then become that 

(a) the limit 

lim -V({p}, /ti(eAi), . . . , n m (e\ m )) (3) 

should be well-defined and non-zero (except possibly in a region of phase-space 
of zero measure). This can be interpreted as a requirement that the soft and 
collinear scaling properties of the observable should be the same regardless of 
whether there is just one, or many emissions. 
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(b) the following two limits should be identical, 



lim lim -V({p}, «i(eAi), . . . , K m (eA m ), K m+ i(eA m+ i)) 

= lim -V({p}, Ki(eAi), . . . , K m (e\ m )) , (4) 

i.e. having taken the limit eq. (|3j). the addition of an extra much softer and/or 
more collinear emission should not affect the value of the of the observable. At 
first sight this closely resembles normal IRC safety, but actually differs critically 
because of the order of the limits on the left-hand side of eq. (j3J). 

These conditions (until now never formulated), which should hold regardless of how 
precisely the /«j(A) vanish as A — > 0, allow one to translate a restriction on the 
ensemble of emissions, V({p}, ki, . . . , k m ) < v, into a restriction on each individual 
emission, V({p}, ki) < v (modulo NLL corrections discussed below). This is necessary 
in order to ensure exponentiation of the LL terms. 1 

3. be continuously global [8] - this means that for a single soft emission, the observ- 
able's parametric dependence on the emission's transverse momentum (with respect 
to the nearest leg) should be independent of the emission direction, d^d^kt In V({p}, k) 
and d^dinkt In V({p}, k) = 0. In practice, this is perhaps the most restrictive of the 
conditions. It avoids the need to analyse possibly quite complicated angular bound- 
aries between regions with different transverse-momentum dependences and calculate 
the corresponding non-global logarithms [8]. It implies a\ = a 2 = . . . = a n = a. 

Given the above conditions, one can derive the following NLL master resummation 
formula for the probability E(-u) that the observable's value is less than v [9]: 

In E(v) = -f^C t n(L) + r' e (L) (in d t - b £ In ^ + B t T 

+ £ In f ^ va+bl ^ + l n S (T(L/a)) + ln^drj, . . . , C n r' n ) 



(5) 



where L = ln 1/v, Ci is the colour factor associated with Born leg t (Cf for a quark and 
Ca for a gluon), and Ei is its energy, Bg accounts for hard collinear splittings and is —3/4 
for quarks and — (HC^ — AT R rif) / {12C a) for gluons, lnc^ = lnc^ + J^ n ^ ln and for 
incoming legs, the fe are the appropriate (Born flavour) parton densities. 

The functions re(L) contain all the LL (and some NLL) terms and are defined by 
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dk\ a s (k t ) . ( k t \ a/b \ f® 2 dk?a s (k t )^Q 



,1 ki it \v l l a QJ J Q 2 V ^ kf it k t 



1 An interesting exercise is to verify that the JADE 3-jet resolution parameter in e + e , which is known 
not to exponentiate [7], is indeed rlRC unsafe. 
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where a s runs at two-loop order and is to be taken in the Bremsstrahlung scheme [10]. 
Exponentiation guarantees that the LL terms of 7> are in the class a™L n+1 . 

All remaining terms are relevant only at NLL accuracy: r' g = di,ri\ T(L) is given by 

^/.x f Q2 dkf a s (k t ) 

T(L) = / -rr^ 11 ■ (7) 

The process dependence associated with large-angle soft radiation is contained in S(T(L/a)), 
whose form depends on the number of legs: 



n = 2 : In S(t) = -t ■ 2C F In 



Qqq' 

Q 



n = 3: ]n.S(t) = -t 



QqgQq'g or i i Qqq' 



c Aln ^w + 2CFln 



n^: ln gW = - f EQln^ + ln^ ge "T /2M r' r/2) . 

Q Tr(HM) 

where Q 2 ab = 2p a .p b and q, q' and g denote the (anti)-quarks and gluon. The n = 2, 3 
formulae apply to e + e~, DIS and Drell-Yan production, while a process such as gg — > 
Higgs + g would simply involve different colour factors. The n = 4 formula applies to 
hadronic dijet production (1 and 2 label the incoming legs). The quantities H, M and T 
are the hard, soft and anomalous dimension matrices of [11] (modulo normalisations and 
our explicit extraction of the factor t from T, see [9]). 

Finally, we examine the factor T . Without it, eq. © corresponds essentially to the 
probability of vetoing all (independent) emissions k with V({p}, k) > v. But given some 
ensemble of emissions fci, . . . , k m that individually satisfy \^({p}, ki) < v, the observable 
may be such that one still has V({p},ki, . . . ,k m ) > v. It is then necessary to apply 
a somewhat stronger veto in order to guarantee V({p}, k±, . . . , k m ) < v. This (and the 
converse situation of V({p}, fcj) > v being allowed in the presence of multiple emissions) is 
accounted for by the NLL function JF, 

T(P> PM-Lnl P'ln V({p},k U ...,k m ) \ 

HR ^ V XP Y R ln max mW ,^),..,n{p},U}/ 

where R' = R'g, R'i = Cgr' t The average is carried out over ensembles of emissions 
generated as follows (cf. section 2 of [12]): first one specifies the value of the maximum of 
the V({p}, ki), say t> max (<^C 1). For each event (ensemble), a random number (m, formally 
infinite) of emissions is generated, according to an independent emission pattern uniform 
in \nk t , rj and (ft, such that on average, below t> max , the density per unit \iaV({p}, k) of 
emissions on leg I is R' e . To ensure a result containing only NLL terms, one takes the result 
in the limit f max — ► 0. Full details, including the derivation and a treatment of subtleties 
associated with the running of the coupling and the recoil momenta, {p} (determined anew 
for each set of emitted momenta), are given elsewhere [9, 12]. 
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We note that attempting to evaluate T for an observable that is rlRC unsafe will yield 
a result that is either ill-defined or improperly behaved for R' 0. This can be thought 
of as analogous to the divergence of NLO terms of a fixed order calculation for observables 
that are IRC unsafe. 

Before proceeding, some remarks on the master formula are in order. As can be verified 
in a straightforward way, we note that: 

• eq. (jSJ is independent of the frame in which one determines the dg, because the 
frame-dependence of the dg is cancelled by that of the Eg; 

• to NLL accuracy, eq. (jEJ is also independent of the choice of hard scale Q; 

• hard emissions collinear to each leg £ are accounted for through the factor Bg in 

eq. (jHJ), and, in the case of radiation from an incoming leg, also through the mod- 

i 

ification of the corresponding parton density factorisation scale from /x F to fi F v a+b i . 
Hard collinear contributions depend only on the combination a + bi, which is to be 
related to the fact that in this region, for an emission k with a fixed energy fraction, 
the observable behaves simply as V({p}, k) ~ (k t /Q) a+be ; 

• finally the continuous globalness of the observable ensures that, to NLL accuracy, 
eq. (J5J) is insensitive to the details of the observable's dependence on large-angle soft 
gluons, the only relevant information being that for any large-angle emission the 
observables scales as V({p}, k) ~ {k t /Q) a . 

Given the above elements, one could imagine a procedure whereby the applicability 
conditions and the parameters of eq. (J1J) are established by hand, analytically, with only 
the T being determined numerically. A related approach was presented in [12], though 
instead of using a master formula, we had to analytically carry out a resummation for 
a 'simplified' version of the full observable — new results were obtained there for three 
observables in e + e~ — > 2 jets. This was already a considerable improvement over the 
traditional, entirely manual resummation approach, which requires a painstaking analysis 
of the observable's dependence on arbitrary numbers of emissions followed by involved 
mathematical procedures to obtain a result which quite often cannot even be expressed in 
closed form (see [13] for a tortuous example). 

However the introduction of a master formula makes it possible to implement a funda- 
mentally new approach. Given a subroutine that calculates the observable for an arbitrary 
set of four momenta, a computer program can carry out the entire resummation: it first 
establishes whether the applicability conditions hold true and determines for each leg i 
the parameters and functions of (JTJ), ag, bg, dg and gg((f>). 2 This is achieved by probing the 
observable with randomly chosen test configurations of soft and collinear emissions, taking 
the asymptotic limit with the help of high precision arithmetic (we choose to use Bailey's 
portable multiple-precision package [14]). 

2 In our current implementation, for technical reasons, ag and bg are restricted be multiples of 1/4, but 
the extension to any power of 1/2 is trivial. 
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If any of the applicability conditions fail to hold (e.g. for the Jade 3-jet resolution 
parameter in e + e~, which is not recursively IRC safe and so does not exponentiate [7]), the 
program does not proceed, i.e. a resummed answer is provided only when the correctness 
of the result is guaranteed to NLL accuracy. 

This method allows one to make an attempt at the resummation of an arbitrary ob- 
servable in a fully automated way, accessible even to non-experts. This is to be compared 
to the standard approach, involving a painstaking (and historically sometimes error-prone) 
manual analysis of the observable, requiring a search for integral transformations (in up to 
5 variables [13]!) to reduce it to a factorised form — this form is then used for the actual 
resummation, after which one evaluates the inverse transforms. 

Usually the two approaches give indistinguishable results, though in some instances 
one or the other may be preferred: for some observables (involving cancellations between 
contributions from different emissions), exponentiation is only partial, resulting in (jSJ) being 
accurate only up to some finite value of R' (typically of order 1) — beyond this point T 
diverges [12] and only the use of the appropriate integral transform method can give a full 
answer (e.g. Drell-Yan p t resummations with a Fourier transform to impact parameter). 
For certain other observables however, the 'factorising' integral transform has yet to be 
found (e.g. the Durham 3-jet resolution parameter) and a numerical approach represents 
the only way of obtaining a resummed answer. 



3 Resummation in hadronic dijets events 

We have verified that our approach reproduces the analytically known results in e + e~ 
and DIS (e.g. [5,8,13]). Here, to demonstrate its feasibility more generally, we show the 
first resummed result for an event shape in hadronic dijet production. Rapid progress 
is currently being made on measurements [15] and fixed-order predictions [16] for such 
observables, with the results showing a clear need for resummations. We shall examine the 
(global) transverse thrust (as opposed to D0's discontinuously global variant [15]), defined 
as: 

T ± = max£#^, (9) 



where the sum runs over all particles in the final state, p± is the momentum transverse to 
the beam direction (rather than to a given leg, denoted by p t ) and n± is the unit transverse 
vector that maximises the projection. 

The transverse thrust has a couple of features worth commenting: firstly, it receives 
non-negligible contributions from emissions nearly collinear to the beams — thus it will be 
sensitive to radiation from the 'underlying event', making it useful for quantitative studies 
of non-perturbative effects that are qualitatively new compared to those examined up to 
now in e + e~ and DIS. Various other observables will be proposed in forthcoming work [9], 
a number of which will be less sensitive to radiation from the incoming legs, providing a 
good degree of complementarity. Secondly, whereas © sums over all particles, experiments 
can only measure up to some maximum rapidity ?7 max . In the presence of such a restriction 
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leg I 


at 


be 


9t{4>) 


de 


(ln^(0)) 


1 


1 





tabulated 


1.02062 


-1.85939 


2 


1 





tabulated 


1.02062 


-1.85939 


3 


1 


1 


sin^ 


1.04167 


-21n(2) 


4 


1 


1 


sin^ 


1.04167 


-21n(2) 



Table 1: Automatically determined leg parameters for t± in hadronic dijet production (in 
a c.o.m. frame with outgoing legs at an angle cos 8 = 0.2). 



it can be shown that the resummation still remains valid for values of v > e -( a + b min)vm^ ^ 
where 6 min is the smaller of the two incoming leg bi values [9] . 

Let us now examine the automated resummation itself: the quantity to be resummed 
is actually t± = 1 — T±, since it is this that vanishes in the Born limit. The observable 
passes all the (automated) applicability tests and table [TJ generated automatically, shows 
the leg properties for a particular reference Born configuration. The different be values for 
incoming and outgoing legs imply different leading logarithmic structures. The azimuthal 
dependence ge{4>) is tabulated and integrated numerically, except in the case of certain 
easily recognisable analytical functions. The function T has a simple analytical form, 3 
however to demonstrate the feasibility of our whole approach we shall show results based 
on a numerically determined T. 

One further step is needed before presenting actual distributions: our master formula 
applies to individual Born configurations. For example in the case of pp — > 2 jets, the 
analysis is carried out with a fixed rapidity for the pair of jets and fixed values of the 
Mandelstam invariants of the underlying hard process. In contrast experimental measure- 
ments integrate over a range of Born configurations. A priori there is no reason for the 
leg parameters or JF to be independent of the configuration and it could be necessary to 
repeat the analysis for a range of configurations. However for most observables, modulo 
certain permutations of momenta (as can be verified automatically), it is only the de that 
depend on the configuration and they are easily redetermined as one integrates over Born 
configurations. 

The resulting distribution for Tj_ is shown in fig. d decomposed into the most relevant 
underlying hard subprocesses, for the Tevatron run II regime (y/s = 1.96TeV). We select 
events containing two outgoing jets with Ej_ > 50GeV and \rj\ < 1.0 and use the CTEQ6M 
parton density set [17], corresponding to a s {Mz) = 0.118. We have set Q — jip — /ir to 
be the Born partonic c.o.m. energy, though in future work we intend to explore a range 
of alternative scales. As is to be expected, channels with lower overall colour charge have 
broader distributions. We note that the different shapes of the various channels constitutes 
information that might be exploitable in fits of parton distributions. Of course detailed 
phenomenological analyses, both for perturbative and non-perturbative quantities, will 

3 It can be automatically established that t± is additive, ki,..., k m ) = J^i V({p}> &»')> implying 

F = e-~f ER '/T(l + R') [5]. 



7 



Ln(T ± ) 



Figure 1: NLL resummed differential distribution D(t±) = dT l (T±)/d\nr± for different 
underlying hard subprocesses. 

also require matching to fixed-order predictions, another step that we leave to future work. 
Here we just remark that resummed results obtained from the master formula are in semi- 
analytical form (fully analytic but for the pure NLL function J 7 ) , so that they can be easily 
expanded to give the fixed-order coefficients needed when matching. 

4 Conclusions 

In this letter we have provided the elements needed for a novel, automated approach to 
general NLL resummation, specifically for the case of continuously global, exponentiable 
(r*fl)-jet final-state observables in the n-jet limit. Results are obtained simply by specifying 
the Born process (and the number of hard partons) and providing the definition of the ob- 
servable to be resummed in the form of a computer routine, similar to the long-established 
practice for fixed-order calculations, and in contrast to the tedious manual approach that 
has been used up to now for resummations. The results are provided in semi-analytical 
form, making it straightforward to obtain the expansions needed for procedures such as 
matching to fixed-order predictions. 

We have demonstrated that the approach can be implemented in practice, by presenting 
automatically generated predictions for the transverse thrust in hadronic dijet production, 
the first event shape to be resummed in this important process. Only concerns for brevity 
prevent us from showing results for a range of other observables and processes, including 
several new observables in hadronic dijet production and jet rates in e + e~ and DIS. 

An open question is whether such an approach, based on the analysis of classes of 
observables can be applied in other resummations contexts, or in the search for higher 
resummation accuracies. We enthusiastically advocate investigations in this direction. 
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